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Abstract. We show that the 1996 outburst of the X- 
ray binary transient system GRO J1655-40 can be ex- 
plained by the standard dwarf-nova type disc instability, 
followed by an episode of enhanced mass transfer from 
the secondary if the mass transfer rate in GRO J1655-40 
is within a factor < 10 of the stability limit. We argue 
that irradiation of the secondary during the onset of the 
outburst driven by the thermal instability in the outer 
disc can increase the mass transfer rate above the min- 
imum value required for stable accretion. This will then 
produce the period of near-constant X-ray emission seen 
in this system. This scenario can also explain the observed 
anti-correlation between the optical and X-ray fluxes. It 
is generally accepted that optical emission in low-mass X- 
ray binaries is produced by irradiation of the outer disc 
by X-rays. There is also strong circumstantial evidence 
that in order for the outer disc to see the irradiating flux, 
it must be warped. Depending on the warp propagation 
mechanism, either a burst of mass from the secondary or 
viscous decay are likely to decrease the degree of warping, 
thereby causing the decrease in the observed optical flux 
while the X-ray flux remains constant or even increases, 
exactly as observed in GRO J1655-40. Finally, the de- 
crease of the disc warping and, therefore, irradiation will 
cause the disc to become unstable once again, terminating 
the outburst. 
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1. Introduction 

The X-ray source GRO J1655-40 has several unusual 
characteristics. Though, like all soft X-ray transients, it 
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is a low-mass binary system, its I.T-S.SMq secondary or- 
biting a 5.5-7.9Mo black hole (Shahbaz et al. 1999) is 
considerably more massive than is typical in such systems 
(see e.g. a review by Tanaka & Shibazaki 1996). Its mass 
and spectral type (F3-F4) imply that the donor star is 
near (or just beyond) the end of its main-sequence life- 
time. In fact Kolb et al. (1997) and Kolb (1998) assert 
that it appears to be crossing the Hertzsprung gap, i.e. 
expanding towards the giant branch, which would imply 
a very large mass transfer rate, Mt ^ 10"'^^ g (King & 
Kolb 1999). Even if the secondary star in GRO J1655-40 
is still on the main sequence, as argued by Regos, Tout & 
Wickramasinghe (1998), the inferred mass transfer rate, 
Mt ^ 10^"^ g s~^, is orders of magnitude greater than Mt 
estimates in other X-ray transients (see e.g. Chen, Shrader 
& Livio 1997). Moreover, the values quoted above are close 
to the critical rate, above which the irradiated accretion 
disc in GRO J1655-40 would be stable with respect to the 
dwarf-nova type instability (van Paradijs 1996; see Dubus 
et al. 1999b for a most recent discussion of the critical 
rate for this instability), that is generally thought to be 
responsible for the transient behaviour of low-mass X-ray 
binary systems. 

The outbursts themselves are also rather atypical. 
Like other Black Hole (Low Mass) X-ray Transient sys- 
tems (BHXTs), GRO J1655-40 has been quiescent for 
more than 30 years before entering the active phase in 
1994. However, contrary to the behaviour of 'conventional' 
BHXTs, the first outburst was followed in 1995 by two 
others which displayed hard X-ray spectra. Finally, af- 
ter the last burst of emission in July/August 1995 the 
system settled into X-ray quiescence, with luminosity of 

« 2 X 10^2 erg 8"^ 

This quiescent state ended around April 25, 1996 when 
a new, soft X-ray outburst began, which is the main sub- 
ject of this paper. Orosz et al. (1997) obtained BVRI pho- 
tometry close to the onset of the outburst and found the 
rise in X-ray flux delayed by 6 days with respect to the op- 
tical increase. Consequently, Hameury et al. (1997) showed 
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that the rise to outburst, and in particular the X-ray 'de- 
lay', are very well described by the dwarf-nova type disc 
instability model (DIM) if the accretion disc is truncated 
at - 5 X 10^ Rs {Rs = 2GM/c^ is the Schwarzschild 
radius). However, the success of this model is put into 
perspective by the subsequent behaviour of the system 
during the outburst, when, after the first local maximum, 
the soft X-ray light curve rose to a higher luminosity 'flar- 
ing' plateau. The present theoretical understanding tells 
us that this type of light-curve cannot be produced based 
on the standard DIM (Hameury et al. 1998). 

The behaviour of the optical luminosity during the 
1996 outburst of GRO J1655-40 also seems to defy the 
generally accepted idea (based on a solid observational 
background) that in Low Mass X-ray Binaries (LMXBs) 
the optical emission from the accretion disc is due to re- 
processing of X-rays in the outer disc. The fact that in 
GRO J1655--40 the optical flux decreases while the X- 
ray flux increases and then fluctuates around an approx- 
imately constant value is difficult to reconcile with the 
X-ray reprocessing model (Hynes et al. 1998). Finally, the 
list of unusual properties of GRO J1655-40 should be com- 
pleted by the presence of a superluminal 'jet' seen during 
the 1994 outburst (Hjellming & Rupen 1995). 

In this paper we argue that the puzzling behaviour of 
this system can be readily reconciled with the DIM, if we 
take into account the fact that the mass transfer rate in 
GRO J1655-40 is likely to be rather close to the stabil- 
ity limit. In this case, irradiation of the secondary during 
an outburst may be enough to increase the mass transfer 
rate by a factor of a few and thereby push the system into 
the stable regime. This scenario and its expected effect on 
the observed X-ray light curve are discussed in detail in 
§11 Furthermore, in §|| we speculate that the reason the 
optical flux does not rise with increasing X-ray luminosity 
is a reduction in the amplitude of the outer disc warp- 
ing. Using a flared planar disc 'approximation' to describe 
a warped disc, we show that as a result, the intercepted 
X-ray flux and therefore the observed optical emission is 
reduced. We further suggest in §^ that the flattening of 
the outer disc and a consequent decrease in X-ray irradi- 
ation may also be responsible for the ultimate end of the 
stable accretion phase, by raising the stability criterion. 
We conclude with a final discussion and a summary in 

2. The X ray Outburst 

Following Shahbaz et al. (1999), in our estimates below we 
assume that GRO J1655-40 contains a TMq black hole 
and that the mass ratio between the primary mass Mi 
and the secondary mass M2 is q = M2/A/1 — 0.33. We 
also adopt a distance to the system of 3.2 kpc (Hjellming 
& Rupen 1995). 

Figure |^(a) shows the progress of the 1996 X-ray out- 
burst of GRO J1655-40. The ASM (2-12 keV) 'soft' X- 
ray flux began to rise 6 days after the start of the op- 



tical outburst and 15 days later attained a maximum 
corresponding to ~ 0.12LEdd (hereafter we use spec- 
tral fits by Sobczak et al. 1999 to convert ASM fluxes 
to bolometric luminosity values), where LEdd — 1-25 x 
lO^^Af/M© erg s~^ is the Eddington luminosity. After 
~ 12 days of a roughly exponential decline which followed 
the first maximum, the soft X-ray flux fell by about 30%. 
It then began to rise again, and after about two months 
from the onset of the outburst, reached a strongly flaring 
plateau with total luminosity varying around ^ O.lTLEdd- 
This flaring state, which Sobczak et al. (1999) identify 
with the very high spectral state, continued roughly until 
day 200 of the X-ray outburst. Over the next ^ 70 days 
the luminosity fell to less than a third of its peak value. 
Note that this local minimum was followed by a ^ 150 
day long reflare, which we do not try to address in this 
paper (but see §p. Our goal here is to apply the (mod- 
ified) DIM to the main 1996 outburst of GRO J1655-40 
since, as discussed above, this model was successfully used 
to describe the beginning of this event. Various type of 're- 
flares' (observed also in dwarf novae) do not have a clear 
explanation in the context of the DIM (see e.g. Hameury, 
Lasota & Warner 2000). One can speculate that the long 
rise-time of the GRO J1655-40's reflare could be due to an 
inside-out outburst from a non-truncated disc, while its 
flat-topped lightcurve suggests a quasi-steady accretion 
phase. However, in view of incomplete optical coverage 
and uncertainties in models discussed below, we feel that 
an attempt at serious modeling of the reflare would be 
pure guesswork. 

It is very suggestive that the e-folding rise time of ^ 1 
day and the e-folding decay time of ~ 35 days, observed in 
the first 27 days of the X-ray outburst of GRO J1655-40, 
are characteristic of the so-called FRED-type light curves 
seen in many X-ray transients (see e.g. Chen, Shrader &; 
Livio 1997). The rise and decay (see Figs. H[b] and |][c]) 
of the optical and UV luminosity also follow the pattern 
of a FRED outburst with an e-folding decline time close 
to ~ 73 days, considerably longer than that for X-rays. 
This behaviour naturally follows from the DIM in which 
the outer regions of the accretion disc are irradiated by 
the X-rays from the inner region (King & Ritter 1998; 
Dubus, Lasota & Hameury 1999a). One could conclude, 
therefore, that GRO J1655-40 was producing a 'normal' 
X-ray transient outburst when, on day 27 a new event 
occurred which stopped the decline of the X-ray flux. In 
what follows, we will argue that the change in the out- 
burst pattern was due to an increase in the mass transfer 
rate, which moved the system parameters into the range 
corresponding to stable accretion for an X-ray irradiated 
disc. 

The minimum mass transfer rate necessary for a stable, 
steady-state accretion in an X-ray irradiated disc may be 
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Fig. 1. X-ray (a), UV (b) and optical (c) light-curves 
(Sobczak et al. 1999; Hynes et al. 1998) of the 1996 out- 
burst of GRO J1655-40. A comparison of the three panels 
clearly shows anti-correlation between the observed opti- 
cal (and UV) fluxes and soft X-ray count rates. The data 
on panels (b) and (c) has not been corrected for interstel- 
lar reddening. Note that the flux scales in panels (b) and 
(c) supersede those on Fig. 3 of Hynes et al. (1998) which 
were in error. 



written as (Dubus et al. 1999b) 
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Here the phenomenological 'parameter' (in reality a func- 
tion of radius and time, see §^) C provides a simple de- 
scription of the disc irradiation properties through: 



4 _ Mc' 



(2) 



where Tjrr is the irradiation temperature. C — 5 x 10^^ 
is the value that represents best the average irradiation 
flux intercepted by the the outer disc regions in LMXBs 
(see Dubus et al. 1999b|]). We assumed that in outburst 



^ Dubus et al. write that this value is used 'by comparison 
with a formula extensively used in the literature' but in fact 
the numerical value of C is deduced from applications of this 
formula to observations. 



the disc outer radius i?out expands to 90% of the primary 
Roche lobe radius (Smak 1999b), so that 



i?ont = 5.3 X 10 



11 f Ml 



\7Mr. 
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cm 
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Let us note, however, that if the outward propagating heat 
front (Hameury et al. 1998) did not reach the outermost 
disc regions, the disc expansion would be negligible, and 
its outer edge would extend only to ~ 70% of the Roche 
lobe radius. 

The critical accretion rate given by Eq. (^ is very un- 
certain, as discussed in Dubus et al. (1999a), and should 
be considered an order of magnitude estimate only. For 
the accretion in an irradiated disc to be stable the mass- 
transfer rate from the secondary has, therefore, to be 
larger than ^ 10^^ g s'^. For GRO J1655-40 this cor- 
responds to Mt ^ O.lA/sdd (AfEdd = -^Edd/O.lc-^). 

Since L ^ 0.2iEdd during the plateau phase of the 
outburst, the observations are entirely consistent with a 
stable disc accreting at M = Mt ~ 0.2MEdd- For our ar- 
guments to work, however, the mass transfer rate from the 
secondary must be smaller than critical value of O.lAfEdd 
quoted above, since otherwise GRO J1655-40 would not 
be a transient, at least according to the DIM. So what is 
known about Mr in this system? 

Estimates of the mass transfer rate in GRO J1655-40 
quoted in the literature differ by three orders of magni- 
tude. At the high end, Kolb et al. (1997) and King & Kolb 
(1999) argue that the secondary star in GRO J 1655-40 is 
crossing the Hertzsprung gap between the main sequence 
and the red giant branch, and is therefore undergoing a 
rapid envelope expansion (on a time scale ~ 10^ years). 
They conclude that this implies a mass transfer rate of 
order M2/(10''y) ^ lO^^g s~^. Obviously, taken at face 
value, this result suggests that GRO J1655-40 cannot be 
a transient system. Regos et al. (1998) show, however, that 
the type F3-F4 donor can still lie on the main sequence 
if convective overshooting at the core-envelope bound- 
ary is taken into account in stellar modeling. The pres- 
ence of such (or equivalent) mixing is expected for main- 
sequence stars in the relevant mass range (see Regos et 
al. 1998 and references therein). Based on this argument, 
Regos et al. (1998) derive a considerably lower value for 
the mass transfer rate in GRO J1655-40, estimating that 
Mt ~ lO^^-lO^Sg s-i « O.l-l.OxM^Jfj. If this estimate 
is correct, the quiescent accretion disc in GRO J1655- 
40 is still unstable, though with the mass transfer rate 
rather close to the critical one. In this case, an increase in 
Mt by a factor of a few during the outburst would stabi- 
lize the accretion, and produce the 'plateau' in the light 
curve, just as observed during the 1996 outburst. Inter- 
estingly, Hameury et al. (1997) showed that the spectrum 
of GRO J 1655-40 in quiescence as well as its rise to out- 
burst in 1996 can be very well modeled using the stan- 
dard DIM operating in a truncated disc, assuming that 
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the mass transfer rate is 2 x lO^^g s~^, within the range 
considered plausible by Regos et al. (1998). 

These relatively high values of Mr are somewhat difh- 
cult to reconcile with the estimate of van Paradijs (1996), 
who used the standard method of deducing the mass trans- 
fer rates in X-ray transients by dividing the mass accreted 
during the outburst by the recurrence time. By estimating 
the total emission observed during the 1994-1995 period 
of activity and assuming a recurrence time of > 30 years 
van Paradijs (1996) obtained Mt ^ 8x lO^^g s~^ (see also 
Menou et al. 1999). However, the validity of this method 
for the recent outburst cycle of GRO J1655-40 is ques- 
tionable, since the observed outbursts clearly do not oc- 
cur with even rough regularity. Between August 1995 and 
April 1996 GRO J1655-40 was in true quiescence so the 
recurrence time for the 1996 outburst was only 9 months. 
The relevant mass transfer rate may then be estimated 
by dividing the mass accreted during the 1996 outburst 
(which was roughly 4 x lO^^'g, assuming the standard 10% 
radiative efficiency) by this recurrence time. The result, 
Mt ~ 2 X lO^^g s^^, is in reasonable agreement with con- 
vective overshooting models of Regos et al. (1998). This, 
of course, does not resolve the problem of the pre-1994 
phase, namely that no outbursts were observed during the 
preceding 30 years. 

Overall we conclude that there is significant evidence 
pointing to the mass transfer rate in GRO J1655-40 be- 
ing rather close to the stability limit of an irradiated disc, 
at least during the period directly preceding the 1996 out- 
burst. Furthermore, we can confidently state that the tran- 
sient nature of this system as well as the shape of the X-ray 
light curve seen in 1996 can be attributed to the dwarf- 
nova type disc instability only if A/t in quiescence was 
below, but not too far from the critical mass transfer rate. 
The former condition is necessary to have outbursts at all, 
and the latter is desirable so that a reasonable increase in 
Mt during the first part of the outburst is sufficient to 
stabilize the disc. In addition, the same value of Mt is 
required to explain all the rise-to-outburst properties. 

An obvious reason for an increase in Mt is X-ray irra- 
diation of the secondary star. According to Phillips, Shah- 
baz & Podsiadlowski (1999), X-ray heating strongly affects 
the vertical structure of the irradiated outer layers of the 
secondary. Assuming an X-ray luminosity of 1.4 x 10'^^ 
erg s~^, they obtain for GRO J1655-40 a (maximum) ra- 
tio of the irradiating to the intrinsic stellar flux ~ 6.6. 
They were considering BATSE (20 -200 keV) X-ray fluxes 
observed during the March 1995 event, whereas during a 
typical outburst most of the energy is emitted in softer 
(2-20 keV) X-rays (Sobczak et al. 1999). At the first max- 
imum in the 1996 outburst, the total X-ray luminosity 
was approximately 8 times larger than the value used by 
Phillips et al. (1999), so the flux ratio would instead be 
~ 50 and the expected effect considerably larger. Thus, 
an increase of the mass-transfer rate due to irradiation 
should be expected. In addition, in a recent seminal ar- 



ticle, Smak (1999a) showed that mass-transfer enhance- 
ments play a fundamental role in dwarf-nova outbursts 
(see §^) and therefore they should also be expected in X- 
ray transients. 

As pointed out by Phillips et al. (1999) X-ray irradi- 
ation does not directly affect the vicinity of the Li point 
but heats up matter located higher in altitude. This will 
induce a delay, because the heated matter has to move to 
the Roche nozzle before falling onto the accreting object. 
Such a delay was observed in dwarf novae; when the heat- 
ing of the secondary by UV and optical emission is clearly 
observed, brightening of the hot spot due to the increased 
mass transfer appears a few days after the outburst max- 
imum (Smak 1995). If the increased mass-transfer rate is 
due to irradiation, this delay must be shorter than the disc 
viscous time, in order to keep Mt approximately constant. 
In the opposite case, the light-curve would instead have 
an exponential shape (Hameury, King & Lasota 1988; Au- 
gusteijn, Kuulkers & Shaham 1993; Hameury, Lasota & 
Hure 1997). 



The disc viscous time is 
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using the expression for midplane temperature given 
for the 'hot' Shakura & Sunyaev (1973) solution. At 
GRO J1655-40's outer disc radius this time is longer than 
the total duration of the 1996 activity (^ 500 days), at 
least for values ofQ;w0.1-0.2 assumed in the DIM. 

In our picture, irradiation of the secondary at the time 
of the first maximum in the X-ray light curve, would in- 
crease the mass transfer rate from the companion. To be 
consistent with observations, this increase should affect 
the accretion rate onto the black hole after about two 
weeks, a considerably shorter time period than tvis above. 
In this case, however, the relevant time scale is not that 
given by Eq. (^, but instead the time it takes for a sur- 
face density front, created by a sudden increase in mass 
transfer, to diffuse towards the inner disc regions. In the 
presence of a warp, matter transfered from the compan- 
ion is fed to the disc not at the outer edge but at radius 
~ ^circ/2 (see e.g. Wijers & Pringle 1999), where 



2.13 X 10 
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is the so called circularization radius corresponding to 
the specific angular momentum of the mass transferred 
through the Li Lagrangian point. The time in which this 
surface density excess will reach the black hole is then 
given by (Hameury et al. 1997) 



f^- < f 
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where 6R is the width of the surface density contrast. A 
density excess with SR/Rdrc ~ 0.1 (see e.g. Frank et al. 
1987) would diffuse to the central black hole in ~ 20 days 
in agreement with our scenario. As we pointed out above, 
this number should be considerably longer than the ex- 
pected delay between irradiation of the secondary and the 
increase in mass transfer, as required to produce a flat 
light-curve. Although it is very difflcult to estimate how 
long it takes for the heated matter to reach the Li point 
(see Hameury et al. 1993, for discussion of a related prob- 
lem), a time of ^ few days is reasonable in the sense that 
it would imply subsonic speeds. 

Our explanation of the 'flat top' X-ray light-curve of 
GRO J1655-40 is very similar to the one proposed by King 
& Cannizzo (1998) for light-curves of Z Cam-type dwarf- 
nova systems (see their Fig. 2). These authors do not 
invoke irradiation to explain the increased mass-transfer 
rate, however, but attribute mass transfer variations to a 
starspot. 

Interestingly, the King & Cannizzo (1998) model is 
more successful in our case than it is in explaining Z 
Cam 'standstills'. In Z Cam the luminosity is observed 
to stick at a constant level halfway down from maximum, 
whereas in light curves produced by King & Cannizzo the 
the 'standstill' is at a luminosity higher than the outburst 
maximum, as in GRO J1655-40. King & Cannizzo expect 
this to be the result of the high mass-transfer enhancement 
factor (6) used in their calculation. In reality, however, the 
mass-transfer rate corresponding to stable disc accretion 
is, by construction, always higher than the maximum out- 
burst mass accretion rate. The reason (see e.g. Hameury 
et al. 1998) is that at outburst maximum, before the cool- 
ing wave begins to propagate, the accretion rate in the 
disc is almost constant, i.e. at the outer disc edge (or at 
the hot disc outer edge) it is close to the accretion rate 
corresponding to the critical surface density (Smin mark- 
ing the end of the 'hot' branch of the S-curve representing 
disc equilibria). A stable mass-transfer rate must be larger 
than this value. 

3. The Optical Light Curve 

There is now a general consensus that optical emission in 
persistent and transient (in outburst) low-mass X-ray bi- 
naries is due to reprocessing of X-rays by the outer disc 
(see e.g. van Paradijs & McClintock 1995). This is also 
the case for GRO J1655-40, since a simple estimate of 
the expected optical flux from a non-irradiated disc which 
reproduces the observed X-ray emission falls short of the 
observed optical flux by roughly an order of magnitude 
(see Fig. ||). Moreover, as we pointed out in the opti- 
cal light-curve is entirely consistent with the decline from 
maximum during a typical FRED-type outburst in an ir- 
radiated disc. Longer UV and optical decline time scales 
(~ 73 days, as compared to ~ 35 days for the X-ray flux) 
are simply due to the fact that as the irradiating flux de- 



clines, the outer disc edge becomes cooler and the peak of 
the emission moves into the optical band, thus compen- 
sating for the decrease in the total emission from the outer 
disc. 

However, this simple picture clearly cannot explain the 
behaviour of the optical flux at times later than 30 
days after the onset of X-ray outburst. At this time. X-ray 
flux begins to increase, while the optical and UV emission 
continues to decline. Here we argue how the scenario for 
the X-ray emission of GRO J1655-40 described in §|| can 
reconcile these observations with the X-ray reprocessing 
origin of the optical emission. 

As shown by Dubus et al. (1999b) the outer regions of a 
planar accretion disc cannot intercept the X-rays emitted 
by a point source located at the midplane. Therefore, in 
order for the outer disc to be irradiated, it must be warped 
(the other possible way for the outer disc to see the X-rays 
- an extended irradiating source - fails to explain why 
the outer disc is effectively geometrically thick while the 
vertical equilibrium implies very thin discs). 

The origin and propagation of warps in accretion discs 
is still an open question (see e.g. Pringle 1999). Here we 
consider two possible regimes: one with low viscosity, when 
the warp propagation relies on sound waves; and another 
with high viscosity, when the warp evolution is driven by 
diffusion. As we shall see, in both cases the warp amplitude 
in GRO J1655-40 is likely to decay during an outburst. 

If the viscosity is low, i.e. if H/R > a, where H is the 
half-thickness of the disc, the warp can propagate as a 
non-dispersive wave at approximately the speed of sound 
(Papaloizou & Lin 1995). This regime can be relevant if 
the outer disc regions are not affected by the propagation 
of the heat front during the outburst. In the standard 
DIM the heat-front passage changes the viscosity param- 
eter a from a low, cold (~ 0.01) to a high, hot (~ 0.1) 
value. Therefore, if the heat front does not reach the outer 
disc it is conceivable (since we don't know the physical 
mechanism supposedly responsible for the change of a) 
that in the outer regions we would still have a ~ 0.01, 
while H/R ^ 0.01 (see e.g. Dubus et al. 1999b). In such 
a case the increased stream of mass transfered from the 
secondary deposits matter moving in the orbital plane at 
~ -Rcirc/2. Since various parts of the disc communicate 
efficiently through sound waves, this new component will 
exert a torque on the outer disc reducing the warp on a 
time scale of a few forced precession periods, where the 
precession period is ~ 40 days (Larwood 1998). 

On the other hand, the whole disc could be in a high 
a state during an outburst. Since in a standard accretion 
disc H/R ^ 0.05 in the outer regions, H/R <C a and 
the warp propagation is driven by viscous processes. The 
relevant viscosity is the one corresponding to the vertical 
shear. The ratio of this kinematic viscosity coefficient to 
the standard (radial) one is approximately 1/2q;^, for a <^ 
1 (Papaloizou & Pringle 1983), and therefore the warp 
damping time is tdamp ~ 4a^ivis- At the outer disc edge of 
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GRO J1655-40, ^damp is then about 100 days for a = 0.1, 
in very good agreement with our scenario. Note that in 
this regime the increase of mass transfer would not affect 
the warp. 

Pringle (1996) found that warp can be radiation 
driven. In such a case the warp's viscous decay could be 
prevented by irradiation. For this to happen, the growth 
rate of the radiative instability must be shorter than the 
viscous damping time tdamp- This condition can be written 
as (e.g. Wijers & Pringle 1999): 



7crit > 3.21 ( -^ 



a 

OA 



M 
R 



1/2 



10 



'1/2 



(7) 



where rj is the accretion efficiency and 7crit ~ 0.1 is the 
critical ratio of the radiative growth to the viscous damp- 
ing times. One can see that for GRO J1655-40 the in- 
equality above can be satisfied only for a values higher 
than the ones usually assumed in the DIM. Unless such 
values are assumed a warp will be viscously damped on a 
time-scale estimated above. 

Whatever the mechanism of warp decay, it would result 
in the reduction of the irradiating flux intercepted by the 
disc, and a consequent decrease in the observed optical 
flux from the system. 

To illustrate this argument we calculated a series of 
optical spectra from uniformly accreting thin discs with 
varying degree of irradiation. The value of the mass ac- 
cretion rate was chosen to reproduce the black body com- 
ponent of X-ray emission (Sobczak et al. 1999). In our 
simple treatment here we specify neither the origin nor 
the structure of the warp. However, since we need some 
description of the photosphere of the disc above the orbital 
plane as a function of radius, we use the prescription. 



{R/Rout) , 



(8) 



chosen by analogy with formulae used in the literature, 
which (despite being based on an incorrect assumption 
about irradiated discs) seem to give a correct empirical 
description of the reprocessed X-ray flux (see Dubus et 
al. 1999b). All other properties of the warped disc, as far 
as irradiation is concerned, are described by C, defined in 
Eq. (^) . We use the prescription for irradiation of the outer 
disc by the inner disc edge (e.g. see Shakura & Sunyaev 
1973; King & Ritter 1998), which combined with Eq. (||) 
above gives 



C = Co 



z/R 



R 

Rout 



4/7 



(9) 



where Cout = C(i?out)- Note that for the disc-disc irra- 
diation geometry, the strength of irradiation is quadratic 
in z/R, since it depends on the projections of both the 
emitting and irradiated annuli. 



Of course the exact dependence of z on the disc ra- 
dius is important in determining the shape of the optical 
spectrum. However, there are many other highly uncer- 
tain quantities in the calculation (e.g. radial profile of the 
albedo in the outer disc, angular distribution of the irradi- 
ating flux, details of radiative transfer in the atmospheres 
of irradiated discs) which all contribute significantly to the 
appearance of the disc in the optical band. In addition, we 
are using a flared planar disc 'approximation' to describe 
a warped disc. Since all these uncertainties are hidden in- 
side our parameter Couti the exact choice of z{R) is not 
very important. As far as we are concerned, Eqs. (^) and 
(^ simply describe the distribution of the irradiating flux 
with radius and do not result from some assumed verti- 
cal disc structure. Note especially that we do not assume 
that the irradiated disc is isothermal or adopt a particular 
value for the disc aspect ratio at the outer edge. 

The resulting spectra computed for different values of 
Cout are shown in figure |^. For comparison we also plotted 
the dereddened spectra of GRO J1655-40 observed (in or- 
der of decreasing flux) on May 14, June 8, June 20, June 
30, and July 22, 1996 with HST/FOS red prism. (Note 
that these spectra correspond to the optical and UV data 
points shown on figure ^) Since the contribution from the 
secondary is quite significant, especially at later times, we 
have subtracted from the data an estimate of the quiescent 
optical spectrum of GRO J1655-40, adjusted according to 
orbital phase (see Hynes et al. 1998 and Hynes 1999 for 
a fuller description of data processing). The broad hump 
centered at logiy ~ 14.8 which remains after subtraction 
in the last two spectra has a similar spectrum to the sec- 
ondary, and we suspect is probably due to a residual con- 
tribution from the companion. This suggests brightening 
of the secondary in outburst due to X-ray heating, sup- 
porting our irradiation scenario. 

Though the model spectra are not intended as a formal 
fit to the data because of many simplifications in our cal- 
culation, the agreement between the five sets of spectra 
is fairly good. Figure || shows qualitatively that by de- 
creasing the degree of warping (as described by Cout), and 
therefore irradiation, of the outer disc, we can mimic the 
observed evolution of GRO J1655-40 in the optical band. 
Note how with decreasing Cout the model optical spectra 
become softer, just as observed, even allowing for some 
residual contribution from the secondary in the data. 

One should keep in mind that all model spectra shown 
in figure ^ were computed for a planar disc with z oc i?^/^ 
and a fixed outer radius, given by Eq. (^). This means 
varying Cout in our model corresponds simply to different 
values of the photospheric height at the outer radius. A 
more realistic description of the effects caused by enhanced 
mass transfer should, of course, include both changes in 
the disc profile (different functional form for z{R)) as well 
as possible changes in the disc outer radius. We feel that 
as long as we use a flared planar disc to approximate the 
effects of the warp, the exact disc shape is beyond the 
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Fig. 2. The observed optical/UV spectra of GRO J 1655- 
40 are plotted together with model spectra from a flared 
disc irradiated by X-rays emitted in the inner parts of the 
disc. The data were corrected for interstellar reddening, 
assuming E{B — V) = 1.2; we also subtracted the contri- 
bution from the secondary. The model spectra were cal- 
culated for M = 0.2MEdd and disc size given by Eq. (^. 
Solid curves correspond to different values of the irradia- 
tion parameter at the outer disc edge, Cout; top to bottom 
Cout = 6 X 10-'', 4 X 10-^3 X 10-^2 x 10-^0.8 x 10"^. 
Dashed curve shows the spectrum of a planar (non- 
irradiated) disc. 

scope of this paper. Similarly, it is difficult to calculate 
the change in the outer radius of the disc from first prin- 
ciples. However, observations and numerical simulations of 
dwarf nova discs (e.g. Smak 1984; 1999b) show that at the 
onset of the outburst i?out expands by about 30% (from 
the canonical 70% to nearly 90% of the Roche radius) for a 
constant mass transfer rate from the companion. By com- 
parison, simulations with enhanced mass transfer (say by 
a factor of 30) simply decrease this value to ^ 25%, so the 
difference is too small to be constrained by the data. 

4. End of the outburst 

Because of many uncertainties in the physics of irradia- 
tion, it is very difficult to infer the time scale of the warp 
decay from the optical data. In both regimes considered 
in the previous section, our estimates of the decay time (a 
few forced precession times in the low viscosity case [Pa- 
paloizou & Terquem 1995; Larwood 1998] or the viscous 
damping time in a high viscosity case) are of the same or- 
der as the duration of the plateau phase of the outburst. 
It is therefore very tempting to speculate that the decay 
of the warp is not only responsible for the decline in the 
observed optical emission, but can also, by reducing heat- 
ing through irradiation, terminate the outburst. Even if 



the mass transfer rate from the secondary remains high, 
a decrease in the irradiating flux can cause an increase in 
the minimum rate required for stable accretion (see Eq. |l]) 
which is sufficient to bring the disc back into the unsta- 
ble regime, and therefore, partially shut off the accretion 
in the inner disc where X-rays are produced. However, 
since the disc will most likely remain close to the stabil- 
ity limit, a very short period of mass accumulation at the 
outer edge is required to trigger the next outburst, just as 
was observed for GRO J1655-40 in 1997 (e.g. Sobczak et 
al. 1999). 

5. Discussion and Conclusions 

We have shown that X-ray and optical observations of 
GRO J1655-40 during its 1996 outburst can be recon- 
ciled with the standard disc instability scenario under the 
assumptions that the outer accretion disc is warped and 
that the primary 'standard' outburst triggers a burst of 
mass transfer from the secondary. In our picture, the en- 
hancement in the mass transfer rate was sufficient to cre- 
ate a period of stable accretion, which was observed as 
a ~ 5 month plateau in the X-ray light curve. We have 
argued that this is entirely plausible, since available ev- 
idence is not inconsistent with the mass transfer rate in 
GRO J1655-40 being close to the stability limit of an ir- 
radiated disc. Thus a rather small enhancement in Mt is 
required to stabilize it. In addition, a burst of mass from 
the secondary is likely to decrease the warp amplitude in 
the outer disc. This effect explains why the optical light 
curve showed a decline as the X-ray flux was increasing, 
and can also account for the spectral evolution of the op- 
tical flux. 

Warped disc physics is rather complicated and calcu- 
lating the thermal-viscous instability in a warped disc 
(not to mention the effects of the enhanced mass trans- 
fer) is out of the question, at present. Because of this, the 
arguments presented in this paper are based on the calcu- 
lations for planar discs and are necessarily of a qualitative 
nature. A lot more work is necessary before we can ven- 
ture to make detailed quantitative comparisons between 
model predictions and what is observed. 

It is also important to emphasize that although the 
overall ratio of X-ray to optical flux (determined by the 
value of C at the outer disc edge) is well constrained by ob- 
servations, the detailed optical spectra of irradiated discs 
depend on many unknown parameters. While the angular 
distribution of irradiating flux may in some cases be de- 
duced from detailed modeling of X-ray spectra, the exact 
shape of the disc as well as the X-ray albedo can only be 
determined from future theoretical calculations. 

Despite these uncertainties, we believe that the pic- 
ture presented here is realistic and is consistent with all 
the available data on the 1996 outburst of GRO J1655-40. 
One should perhaps view our result as a test of the hy- 
pothesis that this outburst was due to a dwarf-nova type 
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disc instability. It is important to note that in onr argu- 
ments we rely on only two major assumptions, which are 
both supported by extensive circumstantial evidence. Disc 
warping emerges as a natural conclusion from studies of 
X-ray irradiated outer discs, and a mass transfer increase 
during the outburst seems equally plausible when we con- 
sider irradiation of the secondary. 

Similar modifications are also required in the appli- 
cation of the DIM to standard (U Gem - type) dwarf 
novae. In these systems there is clear observational evi- 
dence for irradiation-induced enhancement in the rate of 
mass transfer from the secondary (e.g. Smak 1995). Smak 
(1999a) studied the two types of 'inside-out' outbursts 
which are observed in dwarf novae: the 'narrow' and 'wide' 
ones. He showed, that their properties (various widths but 
similar amplitudes) can be explained if during the out- 
burst the mass transfer rate is increased. Narrow outbursts 
correspond to moderate (factor 2) enhancements whereas 
wide outbursts correspond to major enhancements which 
temporarily put the accretion disc into the stable regime. 
It is, therefore, reasonable to expect that the variety of 
light curve types observed in X-ray transient systems may 
be due, at least in part, to mass transfer variations. 

GRO J 1655-40 is an exceptional system but is not the 
only one showing a 'flat top' outburst light curve (see e.g. 
Chen et al. 1997). One may hope, therefore, that such 
systems could be observed in a not too distant future. 
Multi- wavelength observations of such a transient event, 
in particular the correlation between X-ray and optical 
fluxes, could provide a useful test of our scenario. 

At least one feature of the observed behaviour of 
GRO J1655-40 still remains a puzzle. If the mass transfer 
rate in quiescence is always close to the stability limit, it 
is unclear how this system could have sustained a > 30 
years quiescence period. At such a high rate of accretion, 
the mass accumulation period at the outer edge of the disc 
would be very short (on the order of a few months to a 
year), with the successive outbursts following each other 
after a correspondingly short time. In fact, the behaviour 
observed in this system since 1994 should be the norm, 
rather than the exception. One possible explanation can 
be that during some periods accretion in GRO J1655-40 
is non-conservative, e.g. King & Kolb (1999) propose that 
a lot of matter can be lost into a jet. However, more theo- 
retical work as well as observations are required before we 
can assemble a comprehensive picture of this fascinating 
system. 
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